An unconditionally stable algorithm for generalized thermoelasticity based on operator-splitting and time-discontinuous Galerkin finite element methods

An efficient time-stepping algorithm is proposed based on operator-splitting and the space–time discontinuous Galerkin finite element method for problems in the non-classical theory of thermoelasticity. The non-classical theory incorporates three models: the classical theory based on Fourier’s law of heat conduction resulting in a hyperbolic–parabolic coupled system, a non-classical theory of a fully-hyperbolic extension, and a combination of the two. The general problem is split into two contractive sub-problems, namely the mechanical phase and the thermal phase. Each sub-problem is discretized using the space–time discontinuous Galerkin finite element method. The sub-problems are stable which then leads to unconditional stability of the global product algorithm. A number of numerical examples are presented to demonstrate the performance and capability of the method

A p-adaptive, implicit-explicit mixed finite element method for diffusion-reaction problems

A new class of implicit-explicit (IMEX) methods combined with a p-adaptive mixed finite element formulation is proposed to simulate the diffusion of reacting species. Hierarchical polynomial functions are used to construct a conforming base for the flux vectors, and a non-conforming base for the mass concentration of the species. The mixed formulation captures the distinct nonlinearities associated with the flux constitutive equations and the reaction terms. The IMEX method conveniently treats these two sources of nonlinearity implicitly and explicitly, respectively, within a single time-stepping framework. A reliable a posteriori error estimate is proposed and analyzed. A p-adaptive algorithm based on the proposed a posteriori error estimate is also constructed. The combination of the proposed residual-based a posteriori error estimate and hierarchical finite element spaces allows for the formulation of an efficient p-adaptive algorithm. A series of numerical examples demonstrate the performance of the approach for problems involving travelling waves, and possessing discontinuities and singularities. The flexibility of the formulation is illustrated via selected applications in pattern formation and electrophysiology